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Abstract. The B-mode polarization of the cosmic microwave background on large scales 
has been considered as a probe of gravitational waves from the cosmic inflation. Ongoing 
and future experiments will, however, suffer from contamination due to the B-modes of non- 
primordial origins, one of which is the lensing induced B-mode polarization. Subtraction of 
the lensing B-modes, usually referred to as delensing, will be required for further improve¬ 
ment of detection sensitivity of the gravitational waves. In such experiments, knowledge of 
statistical properties of the B-modes after delensing is indispensable to likelihood analysis 
particularly because the lensing B-modes are known to be non-Gaussian. In this paper, we 
study non-Gaussian structure of the delensed B-modes on large scales, comparing it with 
that of the lensing B-modes. In particular, we investigate the power spectrum correlation 
matrix and the probability distribution function (PDE) of the power spectrum amplitude. 
Assuming an experiment in which the quadratic delensing is an almost optimal method, we 
find that delensing reduces correlations of the lensing B-mode power spectra between differ¬ 
ent multipoles, and that the PDE of the power spectrum amplitude is well described as a 
normal distribution function with a variance larger than that in the case of a Gaussian field. 
These features are well captured by an analytic model based on the 4th order Edgeworth 
expansion. As a consequence of the non-Gaussianity, the constraint on the tensor-to-scalar 
ratio after delensing is degraded within approximately a few percent, which depends on the 
multipole range included in the analysis. 
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1 Introduction 

Measurement of the B-mode polarization of the cosmic microwave background (CMB) on 
angular scales larger than a few dozen arcminutes has been considered as the best avenue to 
probe the primordial gravitational waves. The joint analysis of BICEP, Keck and PLANCK 
reported the constraint on the tensor-to-scalar ratio as ro.os < 0.12 (95% C.L.) [1], while the 
much tighter constraint was recently obtained from PLANCK 2015 as ro.oo 2 < 0.08 (95% 

C.L.) [2], where the subscripts 0.05 and 0.002 are the pivot scale of the primordial power 
spectrum in unit of Mpc“^. Still, there is no evidence for the presence of the primordial 
gravitational waves. 

On large scales, Galactic foreground emission, which is one of the most significant con¬ 
taminations, is expected to dominate over the B-mode signals of the primordial gravitational 
waves (sometimes called as the primary B-modes). Many studies have been devoted for 
foreground rejection techniques [3-5] and it would be possible to remove the foreground con¬ 
tamination sufficiently for the detection of the primordial gravitational waves of r ~ 0.001 
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[6]. Even if the foreground contamination is successfully removed, there is still a significant 
contamination due to gravitational leasing from mass distribution between the last scatter¬ 
ing surface and an observer. The gravitational leasing effect on the propagation of CMB 
photons disturbs the spatial pattern of the polarization map, which converts a small portion 
of E-modes into B-modes [7]. The leasing B-modes on large scales has a spectrum like white 
noise and their amplitude is comparable to that of the primary B-modes of r = 0.01 on the 
scales of the recombination bump (£ ~ 10-100) [8]. In addition to precise measurement of the 
leasing B-modes, subtraction of the leasing B-modes, usually referred to as delensing, will 
be required in ongoing and future CMB experiments. It is expected to improve detection 
sensitivity of the primary B-modes [9-13] and even the signals of other non-leasing sources 
such as cosmic strings (e.g. [14, 15]), specific phenomena of modified gravity theories (e.g., 
[16]), and self-ordering scalar fields [17] (more generically, any cosmic defect network [18]). 

For delensing, we need to know the leasing mass distribution, which is described by the 
leasing potential. Measurement of the leasing potential became realized recently by multiple 
observations of CMB polarizations such as ACTPol [19], PLANCK [20], POLARBEAR [21] 
and SPTpol [22, 23]. Near term and next generation CMB observations such as Advanced 
ACT [24], Simons Array [25], SPT-3G [26], CMB Stage-IV [27] will greatly improve sensitivity 
to the leasing potential. In addition, mass tracers at high redshifts such as the cosmic infrared 
background would be also able to be used for delensing [12, 28]. 

To perform likelihood analysis of the B-modes, we should understand and characterize 
statistical properties of the leasing and delensed B-modes. The leasing B-modes are known 
to be non-Gaussian field, and their statistical properties and the impact on the cosmological 
parameter constraints have been discussed in several works [29-32]. These studies showed 
that off-diagonal correlations of the B-mode power spectrum lead to notable degradation 
in the cosmological parameter estimation. On the other hand, statistical properties of the 
delensed B-modes, the knowledge of which is indispensable to analyses in future B-mode 
experiments, has not been explored in detail. 

The main purpose of this paper is to explore how delensing changes the non-Gaussian 
properties of the lensing B-modes, in particular the off-diagonal elements of the power spec¬ 
trum correlation matrix and the probability distribution function (PDF) of the power spec¬ 
trum amplitude. We also discuss its impact on the estimation of the tensor-to-scalar ratio. 
To explore statistical properties of the delensed B-modes, we synthesize Monte Carlo samples 
of simulated CMB maps. Using the simulated maps, we estimate the lensing potentials and 
compute the delensed B-modes. The samples of the lensing and delensed B-modes are used 
for simulating the power spectrum covariance and the PDF of the power spectrum ampli¬ 
tude. To understand the simulation results, treating the delensed B-mode polarization as 
a weakly non-Gaussian field, we construct an analytic model which explains the simulated 
power spectrum covariance and also the simulated PDF of the power spectrum amplitude. 

This paper is organized as follows: In Sec. 2, we describe the procedures of our delensing 
analysis and generation of Monte Carlo samples. In Sec. 3, we show the simulation results, 
focusing on the power spectrum covariance and PDF of power spectrum amplitude of the B- 
mode polarization. In Sec. 4, we present analytic expressions for the statistics of the delensed 
B-modes and compare them with the simulation results. In Sec. 5, we show the impact of 
the non-Gaussianity on the estimation of the tensor-to-scalar ratio. Finally, Sec. 6 is devoted 
to our conclusion and discussion. 

Throughout this paper, we assume a flat ACDM model characterized by six parameters 
which are the baryon density (Db^^); non-relativistic matter density (Dm/i^), dark energy 
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density scalar spectral index (ng), scalar amplitude defined sX k = 0.05Mpc ^ (^s)) 

and reionization optical depth (r). The cosmological parameters have the best-fit values of 
PLANCK 2013 results [33]; = 0.0220, = 0.1409, Ha = 0.6964, Ug = 0.9675, 

Ag = 2.215 X 10-9, and r = 0.0949. 

2 Simulation Method 

In our analysis, to show non-Gaussian structure involved in the delensed B-mode power 
spectrum, we first simulate lensed CMB polarization maps, and then perform quadratic 
leasing reconstruction and delensing. The number of our Monte Carlo samples is 10000. In 
this section, we describe our method of simulating lensed polarization maps, quadratic leasing 
reconstruction and delensing. Note that the simulation method described below corresponds 
to our previous paper [34] . 

2.1 Map simulation 

We denote the polarization anisotropies at a position h on the last scattering surface as 
[Q ± iC/](h). The lensed polarization anisotropies in the direction h, are then given by (e.g., 

[7]): 


[Q ± iU]{ri) = [Q ± iU]{ri + d{n)), (2.1) 

where d is the deflection angle, and is given by the gradient of the leasing potential Viji. 
Here we ignore curl modes in our simulation (see e.g. [35]). Instead of being expressed as 
spin-2 quantities, the following E and B mode polarizations are usually analysed in harmonic 
space (e.g., [7]): 

[E ± iB]em = -Jdn ± 2 Yl^{n)[Q ± iU]{h) , (2.2) 

where we denote the spin-2 spherical harmonics as ±2Yim- Similarly, with the spin-0 spherical 
harmonics, Wm; the leasing potential is transformed into the harmonic space as 

4>lm = (”■)<(>(«.) • (2.3) 

To simulate the lensed CMB polarization map, we first compute the unlensed angular power 
spectra of the E-mode polarization {Cf^) and leasing potential {Cff) with GAME [36]. In 
Lenspix [37], according to these spectra, the harmonic coefficients, Eirn and (piM, are gener¬ 
ated as zero mean random Gaussian fields. Note that, in our leasing simulation, the primary 
B-mode polarization at the last scattering surface is not included since the lensed primary 
B-mode is much smaller than the leasing B-mode. Then, E^m is transformed into Q and U 
maps while 4>im is transformed into </>(«■). Finally, the Q and U maps are remapped according 
to Eq. (2.1) using (p{n). In our simulation, the Healpix pixelization parameter (nside) is set 
to be 2048. 

Expanding Eq. (2.1) up to the first order of the leasing potential, the B modes of the 
lensed polarization field are described as (e.g., [38]) 

= -i E E (L ii m) , (2.4) 

£'m' LM ^ 7 
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where we ignore the primary B-mode. The quantity represents the mode coupling 

induced by the lensing: 


^ +1, + <-«- +1) + i(i +1)] 


fe t o 

V2 -2 oy ■ 

(2.5) 


Here is unity if f + / + L is an even (odd) integer and zero otherwise. Eq. (2.4) 

is known to be a good analytic approximation for the lensing B-modes on large scales. From 
Eq. (2.4), the lensing B-modes in the absence of the primary B-modes are simply expressed 
in terms of a convolution between the E-mode polarization and lensing potential as 

Bim = ^£m[E, c/)] , (2.6) 


where we simplify Eq. (2.4) by defining a convolution operator for two multipole moments: 

^£m[oi, /3] — —i ^ ^ ■ (2-7) 

£'m' LM ^ 7 

2.2 Quadratic lensing reconstrnction 

From the lensed polarization map, we reconstruct the lensing potential </> with the following 
method. Lensing induces the off-diagonal elements of the covariance matrix {(. ^ I' or 
m 7 ^ m') between two lensed polarization anisotropies (X, Y = E,B) as 

= Y.{L!n' m) + 0{cP ‘^), (2.8) 

L,M ^ 7 


where the operation (• • ycMB denotes the ensemble average over the primary CMB anisotropies. 
The weight function in Eq. (2.8) is defined as [39] 


i*EE _ c(~^) /^EE I c(H") /<~rEE /q q\ 

Jee'L — ^ 

=i4ylcr. (2.10) 


Eq. (2.8) motivates the following form for a quadratic lensing estimator (see e.g., [39]): 



* 



1 r jXY 1 * Xem Ye'm' 
A^y \-he'L\ ^xx ^YY 


( 2 . 11 ) 


where X and T are observed polarization anisotropies. The quantities, and Cj^, are 
their angular power spectra. is 2 for XY = EE and 1 for XY = EB. The quantity 
is given by 


A 


XY 

L 


1 lfXYi*fXY 

U££'l 1 h£'L 

2L + l^ 
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Note that, to mitigate biases coming from the higher order terms of (f>, we use the lensed 
power spectrum {Cf^) in Eq. (2.10) [40, 41] instead of the unlensed power spectrum. The 







optimal estimator of the lensing potential is then obtained as a linear combination of the EE 
and EB quadratic estimators: 


with 


^LM - > 


\ _ 

L - (^EE)-I + (^EB)-I 


(2.13) 


(2.14) 


In lensing reconstruction, we add a random white noise with a Gaussian beam factor 
whose power spectrum is described as 






£(£+ 1)02 

81n2 


(2.15) 


Here Tcmb = 2.7K is the mean temperature of CMB, 0 is a beam size, and Ap is a noise level 
of polarization measurement. The multipole £ is an integer between 2 and 2000. The fiducial 
values for the experiment assumed in our analysis are Ap = 6^K-arcmin and 0 = 4 arcmin. In 
such a situation, quadratic lensing reconstruction described here is a nearly optimal method 
[42]. We estimate the lensing potential, using the E and B-modes up to £ = 2000. Note that 
these values are similar to those of the Simons Array [25] which is expected to cover a wide 
range of the entire sky. Advanced ACT also has a similar instrumental specification [24] . 


2.3 Quadratic delensing 

Once we obtain the lensing potential using the quadratic estimator, we can estimate the 
lensing B-modes in a similar manner of Eq. (2.6). On large scales, the lensing B-mode is well 
estimated by [10]: 

(2.16) 

where, again, the operator convolves two fields in the same way that it forms the lensing 
B-modes from the E-modes and lensing potential. We denote EJ^ and as the Wiener 
filtered multipoles, i.e., the measured E-mode polarization and lensing potential multiplied 
by the corresponding Wiener filters {W^ and W^) defined as [9, 10] 


/-yEE 

ttaE _ 

' Cf E + > 

(2.17) 

wf = ^ 

Cf^ + Af 

(2.18) 

The residual field of the B-mode polarization after delensing 

is then evaluated as [10] 


(2.19) 


To avoid delensing bias (see e.g. appendix A of Ref. [34]), only the multipoles between 
£ = 301 and 2000 are taken into account in the calculation of lensing reconstruction. Then, we 
evaluate the delensed B-mode polarization only between £ = 2 and 300. Note that, as shown 
in our previous paper [34] , the filtering of large scale multipoles in lensing reconstruction is 
required also for an unbiased estimate of the lensing potential in presence of 1// noises. 
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Figure 1. Correlation coefficients of the angular power spectrum of the lensing (Left) and delensed 
(Right) B-modes between different multipoles defined in Eq. (3.4). The correlation coefficients are 
computed from 10000 realizations of the lensed and delensed QU maps. We apply flat binning. The 
bin width is Ai = 15 and the multipoles between i = \ and 300 are divided into 20 bins. Note that 
the diagonal elements are removed for an illustrative purpose. 


3 Non-Gaussian Signature of Delensed B-mode Polarization 

3.1 Power spectrum covariance 

Significance of non-Gaussian properties of the B-mode polarization can be seen in the power 
spectrum covariance. In the case of a Gaussian field, the power spectra at different multipoles 
are uncorrelated, and the covariance matrix has only diagonal elements. On the other hand, 
a non-Gaussian field such as the lensing B-modes has off-diagonal elements of notable magni¬ 
tude in the covariance matrix as shown in the previous studies [29-32]. Here, we present the 
simulated covariance matrix of the delensed B-mode power spectrum and compare it with 
that of the lensing B-mode power spectrum. Note that, to see non-Gaussian structure of the 
delensed B-modes, we do not include the tensor perturbations at first. The case including 
the tensor perturbations is discussed in Sec. 5. 

To compute the power spectrum covariance, the multipole is binned into 20 bins. The 
binned angular power spectrum at rth multipole bin is defined as 

= ^ Q, (3.1) 

where Nj is the number of multipoles in zth bin, and we choose ii = 15f -|- 1. For the first bin 
(i = 1), we remove t' = 1. The angular power spectrum is computed with the usual method: 


C,= 


2^-b 




(3.2) 


The power spectrum covariance between each multipole bin i and j is then evaluated as 


Gov,, = (GC,) - (G)(a,). 


(3.3) 
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Figure 2. Histogram of the estimated power spectrum amplitude of the lensing (Left) and delensed 
(Right) B-modes. The number of the samples is 10000. The solid line shows the PDF obtained on 
the assumption that the lensing/delensed B-modes obey Gaussian statistics. We show the ratio of 
the standard deviation computed from the simulation (cr^j„) to that in the case of a Gaussian GMB 
field (erg)- The skewness and kurtosis of the histogram are also shown. The fiducial power spectrum 
is evaluated as the sample mean of the simulated power spectra. 


(• • •) is the sample mean of the operand. 

In Fig. 1, we show the correlation coefficients of the lensing (Left) and delensed (Right) 
B-mode power spectrum defined as [32] 


_ Covij 

Hij — 


i/CoviiCov 


jj 


(3.4) 


The values of the off-diagonal elements increase at higher multipoles. The typical values of 
the correlation coefficients of the lensing B-modes are ~ 6 — 8% at £ ~ 200 — 300, which 
is consistent with the previous results obtained in Ref. [32] On the other hand, we can 
confirm that the off-diagonal elements of the correlation matrix of the delensed B-modes are 
definitely smaller than those of the lensing B-modes. These results imply that the delensing 
operation reduced the off-diagonal covariance more efficiently than the diagonal covariance. 


3.2 Probability distribution function of power spectrnm amplitude 

To discuss the non-Gaussian covariance more quantitatively, next we consider the estimation 
of the power spectrum amplitude. In the case of a Gaussian field, the optimal estimator of 
the power spectrum amplitude is given by [43] 

E.(2^ + l) 

where Ai is the amplitude parameter at each multipole 

Ai = ^. (3.6) 

^ Note that Ref. [32] showed the correlation coefficients with 2 times wider bins in which the off-diagonal 
elements are apparently larger than those in Fig. 1. 


(3.5) 
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The angular power spectrum in the denominator (Ci) is the fiducial power spectrum defined 
as the ensemble average of Q. The parameter Ai is a measured amplitude relative to Ci at 
each multipole. In the case of a Gaussian CMB field, the expected variance of A is given by 


A\2 _ ^ 

E.(2^ + l)' 


(3.7) 


Deviation of the simulated variance from (cr^)^ indicates a non-Gaussian signature. 

The histograms of the estimator A (3.5) of the lensing and delensed B-modes are shown 
in Fig. 2. They are synthesized from 10000 realizations. In the estimation of A, we use the 
B-modes between i = 2 and 300. We evaluate the fiducial power spectrum as the sample 
mean of the simulated 10000 power spectra. Consequently, the sample mean of A is unbiased. 
The standard deviation of A (cr^j^) divided by that in the case of a Gaussian field (u^) is 
shown in the figure. 

Before delensing, the standard deviation of the estimator is ~ 22% larger than that in 
the case of a Gaussian CMB field. On the other hand, after delensing, the discrepancy is 
decreased to ~ 15%, which is considered to be an outcome of the decrease of the off-diagonal 
correlations. Note that, even in the case that the lensing/delensed B-modes are assumed 
to be Gaussian, the skewness and kurtosis of the amplitude estimator have non-zero values 
because the amplitude estimator is quadratic in the CMB field. In our case, the skewness and 
kurtosis of the estimator are 0.0057 and 1.3 x 10“^, respectively. The skewness and kurtosis 
are also shown, and they deviate from those in the case of a Gaussian CMB field within only 
< 0.03 for the skewness and < 0.07 for the kurtosis, respectively. 


4 Modeling Statistical Properties of Delensed B-mode 

To understand the statistical properties of the delensed B-modes shown in the previous 
section, we construct an analytic model of the PDF of the lensing and delensed B-modes 
which is consistent with the simulation results. The simulation results suggest that the 
deviation from a Gaussian field is appeared primarily as the small increase of the variance 
(cj"^)^. Motivated by this fact, we try to express the PDF of the lensing and delensed B-modes 
as a perturbed Gaussian distribution. The Edgeworth expansion, which assumes a weakly 
non-Gaussian PDF, is used for such purpose. We derive a PDF for the lensing and delensed 
B-modes which includes corrections up to 4th order of the Edgeworth expansion. 

4.1 Probability distribution function for weakly non-Gaussian fields 

We denote nL{= +1)) independent Gaussian random variables as aim- The joint PDE 

oi a = {aim} is given by 


Pg(a) 


1 

(27r)’^r/2|C|l/2 """P 



(4.1) 


Here, we denote the covariance of a^^s as C = (aa*). The PDE including the corrections 
from the Edgeworth expansion at 4th order is given by 


P{a) = [I + k{a)]Pg{a), 


(4.2) 
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where the 4th-order correction term k{a) is described as (see e.g. [44]) 

1 d d d d 


k{a) = 


iiirii 


24Pg(a) ^ da, 


^iimi ^^i2m2 da^^Yti^ 


-P,ia) 


(4.3) 


24 ^ 


iilTli 


mim2m3m4 I ®£imi ®4»7i3®£4m4 




(4.4) 


The quantity aem = [C is the inverse-variance filtered multipole. We also denote the 

connected part of the four-point correlation as 


mim2?n.3m4 


= (a^i miO,l2m2^i3m^O,i4,m4)c ■ 

We assume that the covariance matrix of a,m is diagonal: 

This simplifies the PDF as 




2£+l 

, ^ exp 


E 


I I 


L m=—l 


(4.5) 

(4.6) 

(4.7) 


Note that the PDF described above is used to derive the optimal trispectrum estimator of 
the lensing potential power spectrum [45, 46]. 

In the above PDF, we ignore the third order term in the Edgeworth expansion which 
is expressed in terms of the B-mode bispectrum because the bispectrum of the B-mode 
polarization has the odd-parity symmetry and is not generated in the standard cosmology. 

Let us discuss the statistical property of the power spectrum amplitude, A. Using the 
PDF P{a) given by Eq. (4.7), we obtain the mean of the estimator as follows: 

{A) = {{A + Ak{a)))^ = l. (4.8) 


Here, the operation (• • •) means averaging based on the whole PDF P{a), while the operation 
(• • •)g means that on the Gaussian PDF Pg(a). The derivation is shown in appendix A.l. On 
the other hand, from the PDF described in Eq. (4.7), the variance of the amplitude estimator 
is given by (see appendix A. 2) 


(^A\2 — /72\ 1 (2^ -I- 1)(2T -b 1) Covi£i 


riL 


££' 


where the power spectrum covariance, Cov^^/, is expressed as 
2CgC£i _ 1 


Covee' = 


2i+l 


+ 


(2^ + l)(2f+ 1) 


E(-i) 


m+m' rpM'i' 


(4.9) 

(4.10) 


The PDF of the power spectrum amplitude A is given as the chi square distribution with a 
correction term from the kurtosis (see appendix A.3 for derivation): 

- (^g^)^ 


- r(n^/2) ^ ^ ^ ^ 


1 + riL- 


4:(a, 


A\2 


riL 

nL + 2 


- 2 I -b 1 


(4.11) 


Once we obtain the power spectrum covariance, we can actually evaluate the variance of the 
power spectrum amplitude and the PDF. 
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4.2 Expression for B-mode power spectrum covariance 

Next, we discuss mathematical expression for the B-mode power spectrum covariance. At 
hrst, we summarize the derivation of the analytic expression for the power spectrum covari¬ 
ance of the lensing B-modes developed by Refs. [29-32]. The lensing B-mode power spectrum 
on large scales {i ^ 300) is well described by a convolution of the E-mode polarization and 
lensing potential power spectra (see e.g., [10]): 




1 


(4.12) 


£'L 


where siA is given in Eq. (2.5), and we define a convolution operator: 


hi2 


(4.13) 


The fluctuations of the E-mode and lensing potential power spectra contribute to the fluctu¬ 
ation of the lensing B-mode power spectrum as follows [32]: 




a^BB 

£' 




EE I 
£' + 


a^BB 

^ f,r<P<P 

l' ^^£' 


(4.14) 


The connected part of the covariance due to this contribution is approximately 

given by 


L 


dCf^ 2(C'P)2 dC™ 
2L + 1 dC™ 


+E 


L 


dCf^ 2{Cff 

qcH> 2L + 1 ■ 


Here we use 


{5C^5C^) 


<^r'Xr<x 
—-— 

2f-h 1 


(4.15) 

(4.16) 


where X is EE or </>(/). Note that the power spectrum covariance also contains a fully con¬ 
nected term which is not expressed like the r.h.s. of Eq. (4.15), but it has only negligible 
contribution [31]. The approximate expression for the power spectrum covariance is then 
described as [32] 


CoyfP = {Cf^C^^) - {CD{C^^) 

2(C™)2 _ r dcf^ 2(cf^)2 dCf^ 2{cff 

2i+i + dcf® 2L + 1 dcf^ ^ dcf 2L + 1 dcf 


(4.17) 


We now apply the above discussion to derive the power spectrum covariance for the 
delensed B-modes. Assuming that the E-modes are nearly cosmic-variance limited, and 
using the similar procedure described above, we obtain the following expression as 


Cov 


BB,d 

££' 


2iq 


BB,d^2 


2£+l 


-^££' + Yj 


2{Cf^f 2{cl^f 

dcf^ 2L + 1 dcf^ ^ 2L + 1 

(4.18) 


The actual derivation requires cumbersome calculation and is shown in appendix B.l. 
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Figure 3. Correlation coefficients of the lensing (Left) and delensed (Right) B-mode power spectrum. 
The correlation coefficients obtained from the Monte Carlo simulation (solid) are compared with those 
of the analytic model (dashed). For an illustrative purpose, the multipoles up to 300 are binned into 
12 bins. 




Figure 4. Theoretical PDF of the lensing (Left) and delensed (Right) B-mode power spectrum 
amplitude compared with the histogram shown in Fig. 2. We show two theoretical models; the 
analytic model of the PDF given by Eq. (4.11) (solid), and the empirical PDF obtained by fitting a 
normal distribution with the histogram (dashed). Note that, the y-axis is in logarithmic scale for the 
clarification of the effect on the variance. 


4.3 Comparison with simulation results 

Now we compare the model calculations with the results of the Monte Carlo simulation. 

The correlation coefficients of the lensing and delensed B-mode power spectra are shown 
in Fig. 3. The simulated correlation coefficients are compared with those of the analytic 
models described in the previous subsection. The analytic models capture the behavior of 
the simulated correlation coefficients. 

The analytic PDFs of the power spectrum amplitude described in Sec. 4.1 are compared 
with the histograms obtained from the Monte Carlo simulation in Fig. 4. Our analytic 
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Figure 5. Histogram of the tensor-to-scalar ratio estimated from the tensor plus leasing (Left) and 
tensor plus delensed (Right) B-modes. The number of the samples is 10000. The input value of the 
tensor-to-scalar ratio is r = 0.05. The solid line shows the PDF in the case of a Gaussian B-mode. 
The fiducial power spectrum is evaluated as the sample mean of the simulated power spectra. 


model of the delensed B-modes well describes the statistics of the simulated samples. The 
simulated standard deviation of the amplitude parameter is cr^ = 0.005357. By use of the 
analytic formula, the standard deviation is estimated as = 0.005264 which is in agreement 
with the simulation result within 1.8%. We also show the normal distributions htted with 
the histograms. Contrary to the PDFs in the case of a Gaussian field (Fig. 2), the normal 
distributions with the htted variances well capture the behaviors of the simulated histograms, 
which also validates the use of the 4th-order Edgeworth expansion. 


5 Impact on Estimation of Tensor-to-Scalar Ratio 


Here we demonstrate the effect of non-Gaussianity from the lensing/delensed B-modes to the 
estimation of the tensor-to-scalar ratio r. At a single multipole i, the tensor-to-scalar ratio 
which maximizes the likelihood function is obtained by, e.g., differentiating Eq. (3) of [6] in 
terms of r: 


n 


/^BB /^L 
~ 

Of 


(5.1) 


where the quantities C\ and Cj are the observed, lensing/delensed and tensor B-mode 
power spectrum, respectively. Here the spectrum Cj is for r = 1. Under the assumption that 
the covariance of the observed B-mode power spectrum is diagonal, the expected variance of 
the above quantity is given by 


2 2 (C}; + rCj\ 

'"^~2e + l[ CJ ) 

This leads to an optimal estimator of r as 



(5.2) 


(5.3) 
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Figure 6. Standard deviations of the tensor-to-scalar ratio estimator divided by that in the case 
of a Gaussian B-mode (o’gjjjj/trp as a function of the minimum multipole included in the estimation 
(^min) for the tensor plus lensing (dashed black) and tensor plus delensed (solid blue) B-modes. The 
input value of the tensor-to-scalar ratio is 0.05. The maximum multipole in the estimation is fixed to 
300. 


where is the expected variance of r: 



Note that, to see the impact of the non-Gaussian property of the B-modes on the estimation of 
r clearly, we ignore the instrumental noise in the observed B-modes to be used for estimating 
r. The effect of the instrumental noise on the estimation of r is discussed later. The estimated 

tensor-to-scalar ratio would be close to that obtained from the maximum likelihood method 
2 

The histograms of the estimated tensor-to-scalar ratio, which are obtained from 10000 
realizations of the tensor plus lensing and tensor plus delensed B-modes, are shown in Fig. 5. 
The input value of the tensor-to-scalar ratio is 0.05. In both cases, the mean of the estimated 
tensor-to-scalar ratio equals to the input value within Monte Carlo error. Also in both 
cases, we find that the standard deviation of the tensor-to-scalar ratio increases only within 
a few percent (~ 2%) compared with that in the case of a Gaussian B-mode. Unlike the 
constraints on the power spectrum amplitude, the simulation results show that the standard 
deviation obtained in the presence of the lensing-induced non-Gaussian B-modes is 

not significantly different from that in the case of a Gaussian B-mode (o'g). This is because 
well-defined estimators of r usually extracts information mainly from large scale modes which 
earn a most part of the signal-to-noise of the tensor B-modes, while the B-modes on smaller 

^ Note that the estimator r is not precisely the same as that obtained from the maximum likelihood method 
if the tensor-to-scalar ratio is close to zero and the fractional variance of r is large. In our analysis here, the 
estimated values of the tensor-to-scalar ratio are much larger than the standard deviation. 
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scales, where non-Gaussian properties of the lensing/delensed B-modes are significant, are 
down-weighted in the estimators. Note that the value of (crp is 0.00075 (0.00074) in 
the case of the leasing B-modes, and 0.00088 (0.00086) in the case of the delensed B-modes, 
respectively. That is, obtained in the case of the delensed B-modes is decreased by 
~ 15% compared to that in the case of the leasing B-modes, while cjg is decreased by 14% 
after delensing. This means that the delensing efficiency is almost unchanged whether or not 
the non-Gaussian structure of the B-mode polarization is taken into account. 

Let us discuss the case if the instrumental noise is added in the leasing and delensed 
B-modes. If a white noise is added, the B-modes at high multipoles behave like a Gaussian 
field, and the results of the estimation would be more close to those in the case of a Gaussian 
field. On the other hand, if the large scale B-modes are filtered out, the ratio of the standard 
deviations would be increased. In practical situations, the large scale B-modes are 

significantly contaminated by Galactic foreground emission or 1/f noises. In Fig. 6, the ratio 
of the simulated standard deviation to that in the case of a Gaussian field (Csin^/crp is shown 
as a function of the minimum multipole included in the estimation (^min)- As expected, 
the ratio increases for larger values of f'min- Also, the figure compares the constraint on the 
tensor-to-scalar ratio in the case of the tensor plus lensing B-modes with that in the case of 
the tensor plus delensed B-modes. Note that, if i m\r, is increased to be larger than 100, the 
discrepancy between the two cases grows further. 

6 Summary and Discussion 

We have explored the non-Gaussian structure of the delensed B-mode polarization. Assuming 
an experiment of moderate sensitivity where the quadratic delensing is an almost optimal 
method, we find that not only the diagonal elements of the power spectrum covariance matrix 
but also its off-diagonal elements are reduced due to the delensing operation. In particular, 
the correlation coefficients between the power spectra at different multipoles become smaller 
than those of the original lensing B-modes. As a result, the constraint on the amplitude of 
the delensed B-mode power spectrum is not significantly degraded by the non-Gaussianity, 
compared to that on the amplitude of the lensing B-mode power spectrum. We show that 
the PDF of the power spectrum amplitude is broadened due to the non-Gaussianity, which 
leads to the increase of the variance, while the skewness and kurtosis of the PDF are basically 
consistent with those in the case of a Gaussian field. These features are well captured by 
the analytic model based on the 4th order Edgeworth expansion. Also, we show that the 
non-Gaussianity degrades the constraint on the tensor-to-scalar ratio after delensing within 
a few percent. 

In this paper, we assumed the noise level of Ap = 6^K-arcmin and the beam size of 
0 = 4 arcmin. In fact, our analytic model (and mathematical interpretation based on the 
model) works well as long as such moderate sensitivity experiments are assumed. For ex¬ 
ample, we checked that difference of between simulations and the model prediction is 
~ 1% if Ap = 9/rK-arcmin or 7.5/rK-arcmin, and ~ 2% if Ap = 3.5/rK-arcmin. On the 
other hand, decreasing the noise level, we found that the analytic calculations based on the 
formulas described in Sec. 4 eventually began to deviate from the simulation results. The 
ratio jo^ increases after initial decrease in the simulation, while that of the analytic model 
decreases monotonically. For example, if the noise level is 1.5/iK-arcmin, the analytic model 
overestimates the constraint on the power spectrum amplitude by ~ 7%. In absence of in¬ 
strumental noise (i.e., cosmic-variance limited up to £ = 2000), the discrepancy reaches to 
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~ 14%. This result implies that the approximations used in Sec. 4 are no longer valid. In 
such high sensitivity experiments, since quadratic delensing is no longer an optimal method, 
iterative delensing proposed in Ref. [42] would be used for actual delensing analysis. Proper¬ 
ties of the B-modes delensed in such analysis should be explored through another simulation 
(and another analytic modeling) based on the iterative method. It would be our future work. 
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A Statistical Properties of Power Spectrum Amplitude 
A.l Mean 

Let us first derive the mean of the power spectrum amplitude A. The mean of the power 
spectrum amplitude is given by 

(A) = (A}, + (M(a)},. (A.l) 

Since (|a^mP)g = C'r; the first term becomes 

/ T\ _ ^ A 

^ Er(2^+1) ■ ^ 

On the other hand, the second term vanishes. This is because, as described in Eq. (4.3), k{a) 
is a linear combination of four derivatives with respect to a^rn- Therefore, for any nth order 
polynomial f{aim) with n < 4, a quantity {f{a£m)k{a))g vanishes by integration by parts. 
The mean of A then becomes (A) = 1. 

A. 2 Variance 

Next we derive the expression for the variance (4.9), and also Eq. (4.10). The variance of the 
amplitude estimator is defined as 

(,,.4)2 ^ (a2> - 1 = 1)(2<' + l)((AcAe) - 1). (A.3) 

Using the variance in the case of a Gaussian field (cr^)^ defined in Eq. (3.7), the above 
equation becomes Eq. (4.9): 

/ A\2 _ {2£ + 1)(2^^ -|- 1) {CjCji) — CjCji 

nL ^ 2 QQ/ 

_ (2^ + l)(2f' + 1 ) CoYie 

nL ^ 2 CiCi' 

Here, we denote the power spectrum covariance as 

C0V£i> = {C£C£l) — C£Cp . 


(A.4) 

(A.5) 
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The power spectrum covariance is rewritten as 


Co\u' — {C£C£')g — CiCp + {C£C£ik{a))g . (A.6) 

The sum of the first two terms is the disconnected part of the power spectrnm covariance: 

Covf,, = {C£C£,)g - C£C£, = ■ (A.7) 

On the other hand, the connected part of the covariance is related to the trispectrum as (see 
equations from Eq. (A.10) to Eq. (A.14) for the derivation of Eq. (A.9) from Eq. (A. 8)) 

Cov^^, = {C£Cek{a))g 


1 

(2£ + l)(2f+ 1) 

1 

(2£ + l)(2f+ 1) 


^ ^ {^imO‘£rn^i'rn'<^£'rn'^i^))g 

mm' 

(A.8) 

\ ^ / -1 \ m-\-m' n-iiH'i' 

/ V \ -*■/ -^m,—m,m',—m' * 

mm' 

(A.9) 


The sum of the disconnected part (A.7) and connected part (A.9) gives Eq. (4.10). 

The remaining part of this snbsection is devoted for derivation of Eq. (A.9) from 
Eq. (A.8). Using Eq. (4.3), we rewrite Eq. (A.8) as 


1 

(2£+l)(2f+ 1) 


mm' 


1 

(2£+l)(2f+ 1) 

1 

(2£+l)(2f+ 1) 



(A.IO) 


where we define 


K — / dfl (X£YnCL£r^(X£'Yn'^£'\ 


d 


d 


d 


d 


dcLh mi dcii2'm2 m4 


Pg{a). 


(A.ll) 


Introducing an operator Cj which only affects on as Cia£^m.i = perform inte¬ 

gration by parts which gives 


K = -J da [6££^5, 

d d 

X 


l'^mmi{o,£m T Ci(l£jji)ci£i^'Cl£i^i + (£?Ti -f-)- £ TYl 


d 


dOj£2m2 dcL£^Yrn 


Pg{a). 


(A.12) 
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Repeating integration by parts, we find 

K = Jda [5u^ 

^mmi Smm2 (c2 + ci)af m' ^t'm' 


X —^g(«) 

aOjf^^jYiz C'®t4m4 

— do. [S££j^djyimi^££2^mm2i^2 “1“ l^^m'm2i.^(.'m,' d~ C3®£'m') 

“1“ ^Ui^rnmi^Us^mm^icj, + Ci)(5^'£2^m'm2 (®£'m' d~ C2Q.£'m') 

“1“ ^Ul^mm\{,Oj£rn C.l^lm)^(.'(.2^iTi'm2^l'iz^m'mziS'^ *^ 2 )] “1“ {^Ol -f-)- £ ?Ti )] 


d 


X 


<904 




m4 


/i,2 j,/c=(l,2,3,4) 


i E 


Smmh ^iii Smirii tj ^m'rrij Sm^rrik {ch + Ci){cj + Ck)Pg{a) 
Smmfi ^ili ^rarrii £j^m^mj^£'£}^^m'm}^{,£^h “1“ 9” * (A.13) 


h,i,j,k={l,2,3,4) 


The summation with respect to h,i,j and k is applied for every permutaion of (1,2,3,4). 
Using the above equation, we obtain 


(2i + 1)(2£' + 1) i 

'■ ' '■ ' irv) -m ^ 


1 


E /_ 1 \ m+m' rpiii^ i' 

V ^ m,—m,m',—m' * 


(A.14) 


(2£ + i)(2f+ 1) 

ram' 

This equation equals to Eq. (A.9). 

A.3 Probability distribution function of power spectrum amplitude 

Here we derive the PDF of the power spectrum amplitude. A, shown in Eq. (4.11). The PDF 
of A is defined as 


P{A) = jd^^a (5d (^A- 


oc / d”^ a 6 d [A — 


Yllm 


ni 


2t' + l 


[1 + A:(a)] 2 exp 


£ I 12 

\(4£'m.'\ 


m'=—£' 


2Cr 




2e'+i 


I £' 


e"^^ / d^^a 6d{A- 


Yl£m 


ni 


[1 + fc(a)] . (A.15) 


Here (5 d is the Dirac delta function. Introducing variables u^m so that a^m = 
the above relation is simplified as 


P{A) = Pg{A) J d”^u 5 d (1 - \u\‘^) [1 + k{u, A)] 
^Pg{A))C{A), 


(A.16) 
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(A.17) 


where the PDF in the case of a Gaussian field has the dependence as 

^ ^riT - T 

Pg{A) oc 

We define the correction factor as follows: 

JC{A) = = f (5d (l — + f 5d (l — \u\‘^) k{u,A) . (A.18) 

Pg{A) J J 

Let us consider a simplified expression for }C{A). For this purpose, we use the following 
formula 


j d^^x ( 5 d (1 - \x\‘^)fi\x\‘^) = Sul j dr r”^ ^ 5 d (1 - r^)/(r^) 

= ^jds - s)f{s) = , (A.19) 

where 5^^ = jYini^jT) is the surface area of the unit sphere in — 1 dimension. 

From the above equation, the first term of Eq. (A.18) is given by S'n^/2. On the other hand, 
the evaluation of the second term requires a bit complicated calculation. Using the explicit 
expression for the kurtosis contribution (4.4), we first rewrite the second term of Eq. (A. 18) 
as 


1 


E 


^ mim2‘m^m4 





u 



mi '^i2m2 ^Ijim/i 


QnLA5(,-^l^5-mi,m2{ 1) ^%3m3^£4m4 

+ ■ (A.20) 


Note that this is a real number for a given set of a^mS. The above equation consists of the 
three terms: the term proportional to A^, A and independent of A. We first consider the 
term containing A^. To simplify this term, we use the following general formula: 


'y ^ Tijkm ^ 1^1 ) XiXjXf^Xyji 

ijkm 


= I d"-^x 5d (1 - lajp) { + 3 ^ TIj 

^ i¥=j 


ZhY^Ti 


IJtJ • 




(A.21) 


Here, the quantity Tijki is assumed to be unchanged by the exchange of two idices, e.g. 
Tijki = Tjiki- We define h as 


h= d^^x (5d (l — \x\^) 


| 2 | |2 
Xi\ \Xj\ , 


(A.22) 


where i ^ j- Also, we use the fact that (see Sec. A.4 for derivation): 

h = \ [hD (1 - |a;p) \xi\^. 


(A.23) 


- 18 - 





/4 is independent of the indices {i and j) appearing in the r.h.s. of the above equations. The 
quantity 1^ is computed as 


+ I h = jd^^x 5 d (1- 

i i^j 


X\ 


\ 2 


= / d"^£C 5 d (l - |®P) ^ 


5, 


Xi 


n-L 

2 


(A.24) 


By use of Eqs. (A.21) and (A.24), the first term of Eq. (A.20) is given by 

1 _ r 


.47?l4 


hra. 


_ ^ nL ^ (lofmpl at ra 'I )c 


+ 2 16 ^ CfC^/ 

tmi'm' 


(A.25) 


Next we consider the second term of Eq. (A.20). Similarly as in the case of the first term, 
we define the following integral: 


/2 = y d^^u <5d (l — |w|^) UiUj = 5ij Jd'^^u (5d (l — Iwl^) 

This quantity satisfies 

Y.l2= [d^^x6B{l 


Ui 


— \x\ 


E 


2 _ 


Xi\ = 


This reduces the second term of Eq. (A.20) to 
-A 


-T E <^d(1 - 


(A.26) 

(A.27) 


_ A (lu^ml \at'r^' 


2\„ |2\ 
m'\ /c 


CiCt 

Imi'ml 

Finally, the last term becomes 

E! '^mimlrnsniA f d^^U (5d (l — 6i-^i^6-mi,m2^£3i4^-m3,m4{ — ^) 


(A.28) 


mi+m3 


8 

lirrii 

_ ^riL {\aim\^\a£‘ 

“ IfT ^ CiCt 

imi'm' 

Combining these results, we find 




JC(A) = % + 5, 


{nL/{nL + 2)A'^-2A + 1) {\aim\‘^\atm'\‘^) 


2 ' 16 
We finally obtain the PDF of A as 


E 


imi'm' 


CiCt 


(A.29) 


(A.30) 


nr nr t 

P{A) oc A—I 1 + 


i'aL)/{n-L + 2)A^ — 2A + 1 


^(2£ + l)(2/ + l) 


Cov' 


w 


It 


C,Ct 


(A.31) 


Note that the normalization is given by (ni/2)”'^/^/r(re/,/2). 
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A.4 Formula of integration on unit hypersurface 

Finally, we derive the relationship between the following two integrals: 


Jij = jdPx 5d (l - \x\^) 
Ji = j d”® <^D (l — |®P) 


| 2 |™ |2 


Xi \X 


J I ’ 


In our derivation, we use 

F{t) = [(Tx 5d (1 - 1*1^) ^ 


(A.32) 

(A.33) 


(A.34) 


where ,tn) are non-negative parameters. This function is considered as a 

generating function of Jij and Jj, since these quantities are given as a derivative of F{t) in 
terms of ti'. 


d^F{t) 


Jij = lim , 

OtiOtj 


Ji = lim 


d"^F{t) 


t^+o dtf 


The generating function (A.34) is rewritten as 

_ JJ JLg-is(l-|a:|2) g-(qx2-|-t2x|-|—-|-i„x2) 

= [ds j dxi exp [-{ti - is); 

■' ^ i=i 

TT \ 


M - I'Sy 


27r 


n 

i=l 


Substituting the above equation into Eq. (A.35), we find 

Since I 4 = Jij, we obtain Eq. (A.24). 


^Jij — Ji ‘ 


(A.35) 


(A.36) 


(A.37) 


B Delensed B-mode Power Spectrum Covariance 
B.l Analytic expression 

Let us discuss analytic expression for the power spectrum covariance of the delensed B-mode 
polarization. The estimator of the delensed B-mode polarization is given by Eq. (2.19). Eor 
a realization of our simulation, the lensing potential is obtained from the simulated CMB 
map according to Eqs. (2.11) and (2.13). The estimated lensing potential is decomposed into 
the lensing potential and the other remaining term (the noise of the lensing reconstruction) 
as 

— (t^tm T • (I^d) 

From Eq. (2.16), the estimator of the lensing B-mode is described as 

^^„,[E,(r]+mm[E,rF'], (B.2) 
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where the Wiener filtered multipoles 4>J^ and nj^ are defined by multiplying the Wiener 
filter wf in Eq. (2.18) to the and respectively. The dominant contribution of the 
lensing B-mode at large scales is also expressed by the convolution of the E-mode polarization 
and lensing potential as described in Eq. (2.6) [10]. On the assumption that the lensing B- 
mode is given by Eq. (2.6), the delensed B-mode polarization (in the absence of the tensor 
perturbations) is described as a sum of two components: 

= ^im\E, (1 - W^)(P] - mm[E, ri^]. (B.3) 


In the above equations, we assume that the Wiener filter for the E-mode {W^) is unity up 
to £ = 2000, i.e., uncertainty in the E-mode polarization is dominated by its cosmic variance 
and the instrumental noise is negligible. This approximation would be valid for ongoing 
and future high-resolution experiments since their polarization sensitivity will be better than 
0(l)/uK-arcmin also with small beam sizes of a few arcminute. The filter function for the 
lensing potential iwf) depends on the lensing power spectrum of the theoretical model which 
is identified through a parameter estimation procedure in actual analysis. We assume that 
the model is correctly identified in our simulation and the filter function is henceforth treated 
as a quantity independent of Cf‘^. 

The delensed B-mode power spectrum is computed from Eq. (B.3). Under the assump¬ 
tion that the lensing potential is statistically independent from the residual reconstruc¬ 
tion noise the delensed B-mode power spectrum becomes [10] 


~ (1 - 1U‘^)C'^'^], 


(B.4) 

(B.5) 


where is the power spectrum of nirai and, from Eq. (B.4) to (B.5), we assume that 
corresponds to given in Eq. (2.14). Note that, strictly speaking, the lensing potential 
and reconstruction noise are correlated because the reconstruction noise has contributions 
from the lensing potential and lensed E-mode. The impact of the correlation between 
and Uirn would be, however, negligible at least for experiments assumed in this paper, since 
the approximate form of Eq. (B.5) is in good agreement with the simulation results. This 
issue is also discussed in our previous work [34] where we assumed the similar experimental 
specification. 

To derive the power spectrum covariance in the delensed case, we need to compute the 
correlation of Eq. (B.4). Denoting the first and second terms of Eq. (B.4) as 

Df = “^[U™, (1 - 1U‘^)2C'^<^] , (B.6) 

= Ee[C^^, , (B.7) 

we write the covariance of the delensed B-mode power spectrum as the sum of the following 
three terms; 


Tuf = {DfDi) , 

(B.8) 

= (DFDP), 

(B.9) 

= {DfD^,) + ii^ £'). 

(B.IO) 
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The first two terms would be evaluated analytically with the analogy used in the case of the 
leasing B-modes [31, 32], i.e.. 


rpA<l> _ ^ 

L 

= E 


/ 


dDi 


dDi 


+ 


dDi 


dDi 


dlnCf dlnCf 


2L + 1 \ d\nCf^ dlnCf^ 

2 [ dPj dPf, . , dCf^ dCf^ \ 

2L + l\d\nCf^dlnCf^ ^ d\nCf dlnCf j 


(B.ll) 


and 


rri4n _ 

■^U' — 


E 


/ 


dD^ 


dD^, 


2L + 1 I dlnCf^ dlnCf^ 


+ 


{Wt?Nf 


dCf^ dC^^ 
dlnCf din Cf 


(B.12) 


Here we assume that the power spectra and have Gaussian covariance. Note 

that, as discussed in Refs. [20, 48, 49], the covariance between and can have 

off-diagonal elements and is shown that these are basically negligible. The remaining term 
would be approximated as 


rp 2 <l> 2 n _ 


. 1 . 


dDi 


dD^, 


2L + l\dlnCf^ dlnCf^ 


]\r4'4' a^BB a^BB ^ 

+ (1 - » 0. (B.13) 


By combining the above three terms, we find 

2 ( d{Df + D^) d{Df, + D^,) 


I rpAn I rp2(l}2n _ \ 

^£i' + ^U' + . 


2L + 1 


dlnCf^ 


dlnCf^ 
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n 2 




dCf^ dC^^ 


d\nCt‘<' dlnCp 


E 


2L + 1 


d{Df + D^) d{Df, + D^,) 


dlnCf^ 


dlnCf^ 


a^BB 

+ {l-wtf ^ 


dlnCl^ dlnC^^ 




BB,d 


dC, 


BB,d 




^ 2L + 11 9 In Cf E 5 In G™ + q Cf^ d InCf I ' 


From the first to second equation, we assume that corresponds to A^, and from the 
second to third equation, we use Eq. (B.5). The covariance of the delensed B-mode power 
spectrum given above equals to the off-diagonal part of Eq. (4.18). Note that we compute the 
derivative of the delensed B-mode power spectrum with respect to Cf'^ via Eq. (B.5). 
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B.2 Numerical computation 

Here we describe our method of computing the power spectrum covariance of the lensing and 
delensed B-mode. The power spectrum covariance is given by 

Cov®? = ^ C'l+ Covf^/ + Covfg, , (B.15) 


where we denote the connected part of the covariance as 



P E 2(C'™)^ dCp 

^ aC™ 2L + 1 ’ 

L ^ ^ 

(B.16) 


P _ V QCt 2{cff dc,, 

^ dcff" 2L + 1 ■ 

(B.17) 

To evaluate the connected part of the covariance, Cov^ and Cov*^, we 
as 

rewrite the derivatives 


qqee ^eL\y i ’ 

(B.18) 


dCe r/^EEi 

(B.19) 

Here we define 

~ 2(. + l ^2^‘^eLL') ^L' , 

(B.20) 


-ihiM - 21+ 1 ■ 

(B.21) 

Note that 

E,IA, B] = Y1 BL+dA \. 

L L 

(B.22) 

For instance, Cov^/ 

is then given by 



P E dCi 2(C'™)^ dCp 

^dCf^ 2L + 1 dCf^ 



r)/^EE\2 

— ^ '=:4‘ ) "0 [(^00] 

L 2L + 1 

(B.23) 

Denoting 

9^/^EE\2 

) ^0 [^*001 
^ ~ 2L + 1 

(B.24) 

we obtain 

Coyfe,=Ee[F^',C‘^‘^]. 

(B.25) 
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Similarly, we find 


Covf^,=Ee[C^^,G^'] 


(B.26) 


where we define 

_ V L / '^E f/^EEi /x? o^\ 

Gl — 2L + I ^ ■ (B.27) 

The summations in Ei, Ef^ and Efj^ are efficiently evaluated by use of the reduced wigner d 

functions as described in Ref. [10] . The power spectrum covariance of the delensed B-modes 

is also evaluated in the above similar manner. 
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